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I. INTRODUCTION 



The direct interaction of laser fields with nuclei has been considered impossible for a long 
time, most studies being focused on indirect reactions in which parts of the electronic shell or 
plasmas act as secondary particles 111 [2, S]. The nucleus-laser interaction matrix elements are 
small [4], the polarization of nuclei is ineffective, as the populations of the hyperfine levels are 
nearly equal even at very low temperatures, and the energies of the photons are typically several 
orders of magnitude away from those of the nuclear transitions. Meanwhile, synchrotron radiation 
evolved into a versatile source of light which enables one to carry out a great variety of nuclear 
spectroscopy experiments (e.g. Ref. and references therein). Crystal monochromators even 
allow the generation of narrow-bandwidth partially coherent light out of synchrotron radiation, 
such that, e.g., quantum beats could be observed in time spectra of nuclear forward scattering [|7D. 

It is still generally believed that higher power, brilliance and temporal and transverse coher- 
ence at low wavelengths, possibly combined with ultrashort pulses, will open unprecedented per- 
spectives for related experiments in nuclear physics. The most promising candidates for such a 
fourth-generation light source are free electron lasers (FEL) such as the European X-Ray Free- 
Electron Laser (XFEL) |8, i the Linac Coherent Light Source (LCLS) at SLAC H, the SPring- 
8 Compact SASE Source (SCSS) in Japan lllin . the BESSY high-gain harmonic generation FEL 
in Berlin [[l2Q or the fourth generation light source 4GLS in Daresbury [|l3ll (an extensive list with 
currently operating FELs and facilities under development can be found in Ref. H). Recently, 
in view of the forthcoming novel light sources, the direct resonant laser-nucleus interaction has 
been investigated lilSll . considering super- intense laser fields both in the visible and in the x-ray 
frequency regime. The authors of Ref. [llSll have shown that present and upcoming high-frequency 
laser facilities, especially when combined with a moderate acceleration of the target nuclei to 
match photon and transition frequency, do allow for resonant or off -resonant nl6\ laser-nucleus 
interaction. This opens the possibility of optical measurements of nuclear properties such as the 
transition frequency and the nuclear transition matrix element, further developing the field of nu- 
clear quantum optics. In particular, the coherence of the laser light expected from new sources 
such as the European XFEL [8, ^, with photon energies envisaged up to 12.4 keV, is the essen- 
tial feature which may allow to attain nuclear coherence or interference phenomena reminiscent 
of atomic quantum optics, such as nuclear Rabi oscillations, photon echoes or more advanced 
quantum optical schemes llTI] in nuclei. 



2 



While in atomic systems typically electric dipole transitions dominate the light-matter interac- 
tion because of both their relevant transition frequencies and their larger dipole moments, in the 
case of nuclei, this issue is less obvious. Due to the symmetrical structure of nuclei, multipo- 
larities of higher order than electric dipole (El) are the rule, and often electronic transitions of 
the atom, such as internal conversion (IC), compete strongly with radiative transitions of excited 
nuclei lilSll . The spectrum of i^l transitions is limited to few low-lying nuclear excited states with 
small reduced transition probabilities, and giant resonances, at energies of several MeV, which are 
not directly accessible nowadays with the laser and are problematic even when considering the ac- 
celeration of the target nucleus. In analogy to atomic quantum optics, until now only laser-driven 
electric dipole transitions in nuclei have been studied [15, 1^, completely disregarding the large 
choice of electric dipole-forbidden transitions. 

We therefore investigate in this paper electric dipole-forbidden transitions of nuclei interacting 
with super-intense laser fields and show that, unlike in atomic quantum optics, these are suitable 
candidates for resonant coherent driving of nuclei. Since in most cases near-resonant intermediate 
states are not available and the present laser intensities suppress non-resonant multiphoton transi- 
tions, our study is devoted to one-photon transitions. We consider laser-driven transitions between 
the ground and first excited state of suitable stable or long-lived isotopes, including all the mag- 
netic sublevels in a multilevel model of the nuclear system interacting with the laser field. In order 
to account for dipole-forbidden transitions, the plane waves describing the electromagnetic field 
are expanded in spherical waves characterized by well-defined multipolarity and parity. The nu- 
clear interaction matrix element is expressed with the help of the reduced transition probabilities, 
for which we use experimental values. The interaction matrix elements and population inversion 
for electric and magnetic dipole and electric quadrupole transitions are compared for several sta- 
ble or long-lived nuclei with low-lying first excited states. This comparison includes both cases 
in which the direct laser-nucleus interaction is by itself possible as well as nuclear systems for 
which moderate acceleration of the target nuclei is required to match nuclear transition and photon 
energies. We find that Ml transitions are prospective candidates for resonant coherent driving 
of nuclei due to their large transition probabilities and suitable energies. A comparison between 
electric dipole- allowed and forbidden transitions shows that, unlike in atomic cases, for similar 
transition energies the interaction matrix elements are often of the same order of magnitude. This 
result considerably expands the choice of suitable isotopes since electric dipole transitions from 
the ground state in the keV region are very rare. 



3 



This paper is organized as follows: in Section II, we derive the laser-nucleus interaction Hamil- 
tonian used in this study and calculate the interaction matrix element for higher multipoles, fo- 
cussing on E2 and Ml transitions. Section III presents the numerical results for the interaction 
matrix elements and population inversion of the nuclear multilevel system and discusses aspects 
of a possible experimental observation of the laser driving of nuclei. We conclude with a short 
summary. 



n. THEORETICAL FORMALISM 



For applications in quantum optics, the density matrix theory has been utilized successfully to 
describe the dynamics of the atomic system interaction with laser fields including relaxation lIlTn . 
In the following, we apply the density matrix formalism to the few-level nuclear system driven by 
the electromagnetic field. We use a semiclassical description in which the nucleus is treated as 
a quantum few-level system and the field is treated classically. The equations of motion for the 
density matrix are derived in Section HTAl The laser-nucleus interaction Hamiltonian is detailed in 
Section IIIBl where the electromagnetic field is expanded in spherical waves taking into account 
the specific terms responsible for the considered dipole-forbidden transitions. 



A. Basic quantum dynamics 

The interaction Hamiltonian between the nucleus and the electromagnetic field can be written 
in the semiclassical approximation as 

Hi = -\j d'rj{f)-A{r,t), (1) 

where j(r) represents the nuclear charge current, c is the speed of light, and the vector potential of 
the field in the Coulomb gauge is given by 

Air, t) = —Eke-'^^'e-'^-'%^ + c.c. (2) 

UJk 

The electric field amplitude is denoted by Ek, and the photons are characterized by the wave 
vector k, the frequency uJk, and the polarization a. Furthermore, e^^ is the polarization vector. As 
an example we consider circularly polarized light, so that cr = ±1. The case of linearly polarized 
light can be treated in a similar manner. 
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We restrict the Hilbert space of the nucleus to the excited state |e) and the ground state \g) 
characterized by the angular momenta Jg and Ig, respectively, including their magnetic sublevels 
Me and Mg. The nuclei assumed to be initially in the state \g) are irradiated with the intensity 
X{t) starting from time t = 0. The dynamics of the density matrix p is determined by the Master 
equation flV"] 

d 

th—p=[Ho + Hi,p]+Cp, (3) 

where Hq denotes the unperturbed nuclear Hamiltonian operator and the Lindblad operator C 
describes the various spontaneous relaxation channels. 

The density matrix element corresponding to the states i,j E {e,g}, that are characterized by 
the magnetic quantum numbers Mj and Mj, respectively, is denoted by pij. The optical Bloch 
equations evaluate to 

^Pgg{Mg) = - |lm lj2pUMg,M,)e'--'{h,M,\Hj\Ig,Mg) 



Me 



+ J2^s{Mg,M,)pUMe 



Me 



g-^Pge{Mg, Me) = - iApg,{Mg, M + {Pgg^Mg) " P ^e^M ,)) C-'^^^^^Ig, Mg\Hl\h, M 

_ 7.(Mg,Me) ^^^^^^^ _ ^^^^p^^^M^^ ^ 

^Pee(Me) = ^Im p^e (M„ Me)e^-^* ( Je, Me | i/, | Jg, M^) 

-pUMe)J2js{Mg,M,), (4) 

where the rapidly oscillating terms in the off-diagonal elements have been eliminated by moving to 
a rotating frame. The spontaneous decay of the excited level |e) depends on the magnetic sublevels 
of the nuclear states. With Tg being the total decay rate of the |e) state, summed up over all the 
possible values of Mg and Mg, the partial decay rate 7s (M^, Me) is given by the Clebsch-Gordan 
coefficients C{ji j2 js; "^2 "^3) via 

7,(M„ Me) = ^|^[C(/<; le L' Mg - M, M)]2 r, , (5) 

where L and M denote the photon multipolarity and total angular momentum projection, respec- 
tively. In the optical Bloch equations (HI), 'j^ec stands for an additional dephasing rate to model 
laser field pulses with limited coherence times, while A = luq — ujk '^s the detuning of the laser 



frequency with respect to the nuclear transition frequency luq. The system of differential equations 
with the given initial conditions can be solved once the matrix element of the interaction Hamilto- 
nian is known. In the following section we calculate the interaction matrix element, which is also 
contained in the expression of the Rabi frequency of the system. 



B. Interaction matrix element 



In atomic physics, the electromagnetic field is usually described by the plane-wave expansion, 
such that the vector potential has the expression given in Eq. Q. The photons are characterized 
by the wave number k, the propagation direction k and the polarization a. In the case of nuclear 
transitions, however, it is more convenient to describe the photons in terms of multipolarity and 
parity. We therefore expand the plane wave in multipoles to account for photons characterized by 
the angular momentum L, its projection M and the parity A. We follow the formalism described 



in [|l9|] and express the product in the expression of the interaction Hamiltonian in Eq. ([U) with the 



help of the electric £ and magnetic Ai multipole fields, 

LM 

where D^j -a(k) is the rotation matrix (characterized by the photon angular momentum projection 
M and polarization a) associated to the rotation of the coordinate system which transforms the 
direction of the 2;-axis to k. The magnetic A^j^j{r) and electric ^ijvf(^ multipole fields are given 
by 



A^m{^ = 3L{kr)Y^,{0.^). 

ALA^ = J^^jL-i{kr)Yli^,{e,v) - J^^jL+i{kr)Y,'l^,{e,v) 



= -^Vx ^n{kr)Y,^[{e,cp)) , (7) 

where iL{kr) are the spherical Bessel functions and Yll{9, y?) denote the vector spherical harmon- 
ics. Typically, only a few terms in the sum in Eq. i.e., one or two multipole orders, are needed 
for the accurate calculation of the interaction matrix element. In the case of magnetic multipole 
transitions, the interaction Hamiltonian has the form 

Hi = -^e-^-'=*v^ V(-^)^v/2m/^i,_(^) [ d\jL{kr)^{r) ■ YlUe,^) . (8) 

" LM 



We consider in the following that the direction of the incoming photons is parallel to that of the 
z-axis, so that D^,j_^(k) = 5m -a, where 5ij is the Kronecker delta symbol. As the wavelength 
of the radiation is large compared to the nuclear radius (kRo <^ 1), we use the long wavelength 
approximation, keeping only the first order term (kr)^ in an expansion in (kr) in the spherical 
Bessel function. The interaction matrix element then becomes 



k 



L-l 



(9) 



where -a is the nuclear magnetic multipole moment operator corresponding to the total angular 



momentum L and its projection —a, defined as [|20|] 



Mlm = / rfV(r-*x j(f)) ■ V(r^FiM(^,</^)) • (10) 



c(L + 1 

The symbol ! ! in Eq. ^ denotes the double factorial given by n! ! = n(?i — 2) (n — 4) ■ ■ ■ where 
Kn is 1 for odd n and 2 for even n. Since the multipole moments are spherical tensor operators, we 
can factor out the magnetic quantum number dependence of their matrix elements with the help of 



the Wigner-Eckart theorem [12 ih . 



{h,M,\ML^M,,Mg) = ^j^=^C{hIgL;M, - M, - a)(/e||Mi||/,) . (11) 

V 2L + 1 

The reduced matrix element of the magnetic multipole moment {Ie\\ML\\Ig) is now independent 
of the angular momentum substructure and can be related to the reduced transition probability B 



B{ML,I, ^ Je) = ^j^mML\\I,)\' . (12) 



The interaction matrix element for a certain multipolarity L can then be written as 



X y/2Ig + l^B{Mmg^Q, (13) 

omitting an overall phase factor which is irrelevant, since the population dynamics of the few-level 
system depends on the squared modulus of the interaction matrix element. In a similar manner, by 
identifying the nuclear electric multipole operator in the expression of the interaction Hamiltonian 
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for electric transitions, we obtain for the matrix element the expression 



2/, + K/5(^L,/g^/e). (14) 



ni. RESULTS AND DISCUSSION 

The interaction matrix elements and population inversion for electric and magnetic dipole and 
electric quadrupole transitions are calculated and compared for several stable or long-lived nuclei 
with low-lying first excited states. This comparison includes both cases in which the direct laser- 
nucleus interaction is by itself possible (the nuclear transition energy E allows for direct resonant 
interaction, E < 12.4 keV), as well as nuclear systems for which moderate acceleration of the 
target nuclei is required to match nuclear transition and photon energies {E > 12.4 keV). Signal 
rates are calculated envisaging realistic experimental parameters of the laser and nuclear targets. 
An overview of the important light source parameters is given in Section HIlAi while the numerical 
results are present and discussed in Section UlTBl 

A. Light sources 

Powerful high-frequency coherent light sources are the key element for achieving direct laser- 
nucleus interactions. In the x-ray energy range, FEL is currently the only suitable lasing mecha- 
nism, exploiting, for example, the self-amplified spontaneous emission (SASE) of free electrons 
accelerated and propagated through an undulator. Since the electrons in the FEL are not bound 
to atoms and thus not limited to specific transitions, the wavelength of the FEL is tunable over a 
wide range depending on accelerator energy and undulator parameters. At the future European 
XFEL, photon energies up to 12.4 keV are envisaged, with an average spectral brightness [[22] 
of 1.6x l_0^^j)hotons/(sec-mrad^-mm^-0.1%bandwidth) in pulses of 100 fs at a repetition rate of 
The SPring-8 Compact SASE Source (SCSS) in Japan aims at a slightly lower 



40 kHz m 
spectra 
SLAC 



brightness, with similar photon energies [llin . The Linac Coherent Light Source (LCLS) at 
I I0I1 will deliver photons with energies up to 8.2 keV with an average spectral brightness of 



2.7x10^^ photons/(sec-mrad^-mm^-0.1%bandwidth). Other FEL light sources such as the BESSY 



high-gain harmonic generation FEL in Berlin il2|] or the 4th generation light source 4GLS in 
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Daresbury 11311 can reach only a lower photon energy, up to 1 keV or 100 eV, respectively. As the 
number of laser-driven nuclear transitions per pulse is small, the pulse repetition rate is a crucial 
value for the magnitude of the nuclear excitation effect to be observed. The European XFEL has 
the largest repetition rate of 40 kHz, while the SCSS at Spring-8 and the LCLS PEL at SLAC have 
only 60 Hz and 120 Hz, respectively. The light sources with lower photon energy at BESSY and 
4GLS can provide pulses with a repetition rate of approximately IkHz. 

Since the European XFEL is designed to reach the highest photon energy, pulse repetition rate 
and average brilliance, we use its parameters to numerically estimate the possible magnitude of 
direct laser-nucleus interaction. The peak intensity of the laser can be obtained from the peak 
power divided by the area of the focal spot, X = 6 x 10^^ W/cm^. For a pulse of 10^^ photons 
with an energy of 12.4 keV, the total energy spread of about 10 eV is given by the spectral 
bandwidth of 0.08 %. The nuclear transition width is many orders of magnitude smaller than the 
photon energy spread, implying that only a fraction of the incoming photons will actually fulfill the 
resonance condition. From the flux of photons resonant within the transition width of the excited 
state Fat, we deduce the effective laser intensity Ics, 

Jeff = ^J. (15) 

As far as the magnitude of the peak intensity is concerned, we distinguish two cases. In the first 
case, we consider an optimization of the number of signal photons scattered off of the nuclei 
for a resonant laser field. Since the intensity of the considered light sources typically does not 
allow for full Rabi flopping in the laser-nuclei interaction due to the small width of the nuclear 
resonances, high average intensities rather than peak intensities are desirable. For laser sources 



strong enough to achieve an inversion of the nuclear population (via a so-called vr-pulse nll\\ ). 
the intensity should be adjusted such that the population evolution does not exceed the complete 
inversion, as otherwise the signal photon rate is reduced. Second, if quantum optical scenarios 
which require a Rabi frequency comparable to the radiative decay rates are considered, high peak 
intensities are mandatory. The feasibility of quantum optical schemes with present laser sources 
will be discussed in more detail in the following section. 

If the target nuclei are accelerated, the boosted electric field and the laser frequency i^at are 
given by 

En = {l+PhEL, 

UN = {l+P)lUL. (16) 



where the subscripts and L stand for the rest frame of the nucleus and the laboratory frame, 
respectively. The reduced velocity [3 and the Lorenz factor 7 are chosen such that the laser and 
nuclear transition frequencies match, e.g. a nuclear transition of 100 keV with 12.4 keV laser 
photons correspond to a Lorentz factor of 7 = 4.1. In the boosted reference frame, the effective 
laser intensity is multiplied by the square of the boost factor, 

:^eff = (l + /5)Vxifr. (17) 

For the particular case of very low-lying first nuclear excited states in the eV regime, such as 
the case of the isomeric state of go^Th which will be discussed later in this section, conventional 



light sources as well as FEL lasers can be used. Among the conventional UV lasers, Vulcan ll23h 



at the Rutherford Appleton Laboratory and PHELIX [|24l 12511 at the GSI in Darmstadt can reach 
intensities of up to X = 10^^ W/cm^. The drawback of these light sources is the low repetition rate 
of order of one pulse per 10 min. A better repetition rate of the order of MHz can be achieved, 

n 

e.g. by the VUV-FEL at the 4GLS in Daresbury [130 with a photon energy of up to 10 eV, and an 
average brightness of 5x10^^ photons/(sec-mrad^-mm^-0.1%bandwidth). 

B. Numerical results 

We calculate the interaction matrix elements for a number of nuclear systems involving E2 and 
Ml transitions between the ground and first excited nuclear states. Particular interest is devoted 
to nuclei in which the energy of the first excited state allows for a direct interaction with the laser 
without requiring the acceleration of the target nucleus for the given European XFEL parameters. 
Namely, we investigate the cases of the sgKr, sy^La, g|^Sm, gg^Tm and yg^Os isotopes, which have 
first excited states lying below 12.4 keV. These isotopes have stable or long-lived ground states, so 
that a reduction of the full level structure to the nuclear ground state and the first excited state in a 
few-level approximation is well justified. Each of these states is split into its respective magnetic 
sublevels. For comparison, a low-lying E\ transition at 6.238 keV in yl^Ta is also investigated. 
For the moment, in the framework of a general discussion, we assume that the spontaneous decay 
of the excited nuclear state occurs only radiatively. This corresponds to nuclei in bare ions or cases 
in which nuclear de-excitation via IC is not possible. The radiative decay rate Vg can be calculated 



using the reduced transition probabilities 1I20I1 . 
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where A stands for electric 8 or magnetic M., and E denotes the transition energy. Values for the 
reduced transition probabilities, the interaction matrix elements and the radiative decay rates are 
given in Table HI The reduced interaction matrix element denotes the quantity 



{Ie\\Hj\\I,) = EeffV^y L (2L + 1)!! V^^\/BiML, I, ^ Q (19) 



for the case of magnetic transitions and 



{Ie\m\I,) = E,W27i^ + ^l^^ ^^^^ y/2I, + 1^B{£L,I, /,) (20) 

for electric transitions. Here, Ecs is the effective electric field amplitude, corresponding to the 
effective laser intensity Xcfj. For these cases, the XFEL is assumed to deliver photons in resonance 
with the nuclear transition. The energy levels and the experimental values for B(\L, 1^ — * Ig) 



were taken from [126 . 



27 



30 



3lll . The reduced transition probabilities for the emission and 



the absorption of a gamma ray, respectively, are related through the formula 

BiXL, Ig ^ Q = '^±^B{XL, h ^ /,) . (21) 

In the cases of the IgKr, g2^Sm, gg^Tm and yg^Os isotopes, we neglect the weak multipole mixing 
of the nuclear transitions, considering only the dominant Ml component. The multipole mixing 
ratio values given in Ref. [32] support this approximation, which affects with less than one percent 
the accuracy of the calculated signal rate. 

The values for the reduced interaction matrix element presented in Table IJ show that for low- 
lying nuclear levels with energies of approximately 10 keV, Ml transitions are well-suited for 
direct interaction with strong laser fields. The largest interaction matrix element occurs in the case 
of the 9.746 keV Ml transition of \fOs. 

In order to extend the comparison between electric dipole allowed and forbidden transitions 
we also consider typical examples of stable or long-lived isotopes with the first excited state at 
energies covering the range from 12.4 keV up to 100 keV. Since the width of the first excited 
state determines the effective intensity of the laser field, we chose for several energy ranges and 
transition types isotopes with short-lived excited levels. The effective intensity of the laser Xcff 
is calculated according to Eq. (flTI) . assuming that the XFEL delivers 12.4 keV photons that are 
boosted in the nuclear reference frame to reach the transition frequency. In Table HI] we present 
reduced interaction matrix elements and spontaneous decay rates for the El transitions of Sm, 
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TABLE I: Reduced interaction matrix elements and spontaneous decay rates for a choice of isotopes with 
first excited states lying below 12.4 keV. In the fifth and sixth columns, E and B(XL, Ig — > I^) denote the 



nuclear excitation energy and the reduced transition probability, respectively, taken from Refs. ||26 , 



27 



28 



30, 



Isotope 




le 


L 


E (keV) 


B{\L,Ig^h) (eW^) 


{h\\Hl\\Ig) 


(eV) 


r, (i/s) 


83t7- 
36 -"^^ 


9/2+ 


7/2- 


M1+E2 


9.396 


1.49 X IQ-"^ 


1.61 X 10' 


-10 


2.51 X 10^ 


137t „ 
57 


7/2+ 


5/2+ 


Ml 


10.56 


2.51 X IQ-s 


2.50 X 10' 




6.37 X 10^ 




5/2~ 


3/2- 


M1+E2 


4.821 


7.74 X IQ-s 


7.07 X 10' 


-11 


2.10 X 10^ 




1/2+ 


3/2+ 


M1+E2 


8.410 


1.24 X 10-3 


3.79 X 10' 


-10 


5.96 X 10^ 


i87r)„ 

76 


1/2- 


3/2- 


Ml(+E2) 


9.746 


1.26 X 10-3 


3.84 X 10' 


-10 


9.43 X 10^ 


73 """Ta 


7/2+ 


9/2- 


El 


6.238 


5.18 X 10-6 


6.84 X 10' 


-12 


1.59 X 10^ 



Eu and Dy, the 0+ 2+ E2 transitions of Gd, Dy and U and finally the Ml tran- 
sitions of 7o^Yb, Ho and gl^Eu. All these isotopes have stable ground states, except for 62^Sm 



which has a mean lifetime of r = 46.27 h 1321]. and the very long-lived g2^U ground state with 
mean-life of approximately 4x10^ yr. The lifetime of the nuclear ground state is important when 
considering possible experimental observation of the coherent driving of nuclei, because it affects 
the sample preparation. While the interaction matrix element is the largest for the El transition in 
g2^Sm, this isotope has the disadvantage of its unstable ground state. On the other hand, the Ml 
transition in the stable gy^Ho isotope can be a good candidate for the experimental observation 
of nuclear quantum optics effects. The higher order term in the wave number k causes the E2 
transitions to have smaller interaction matrix elements, as shown in Table HIl 

By solving the optical Bloch equations in Eq. (H)) we can obtain the population of the excited 
state after one pulse, considering as initial condition that the ground state population is equally 



distributed between the corresponding magnetic substates pSi] . In the scenario of laser interacting 
with bare nuclei, after the laser pulse, the excited nuclei decay back to the ground state with the 
emission of a signal photon. Thus, the sum pee of the excited state population pee{Me) over the 
possible magnetic sublevels Me is equal to the (unpolarized) signal photon rate per nucleus per 
pulse A'^giGNAL for the driven nuclear transition. We present in Figure [H the signal photon rate for 
the isotopes considered in Tables IJ and |II] for nuclear transitions with energies above and below 
12.4 keV using the European XFEL laser parameters. The highest signal rate is iVsicNAL = 
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TABLE II: Reduced interaction matrix elements and spontaneous decay rates for a choice of isotopes with 
low-lying first excited states above 12.4 keV. In the fifth and sixth columns, E and B{XL, Ig — > Ig) denote 



the nuclear excitation energy and the reduced transition probability, respectively, taken from Refs. 133 , 
35 



34 



361, 



Isotope 


y 


h 


L 


E (keV) 


B{XL,Ig^Q (e^fm^^) 




(eV) 


Ts (1/s) 


ef Sm 


3/2+ 


3/2- 


El 


35.843 


> 3.50 X 10-2 


2.27 X 10 


-7 


2.55 X 10^ 


63 


5/2+ 


5/2- 


El 


97.429 


1 80 X 10-'^ 


1 06 X 10 


-7 


2 64 X 10^ 




5/2+ 


5/2- 


El 


25.651 


2.65 X 10"^ 


1.08 X 10 


-9 


7.09 X 10^ 


ifGd 


0+ 


2+ 


E2 


88.966 


4.64 X 10^ 


3.27 X 10 


-9 


6.32 X 10^ 




0+ 


2+ 


E2 


80.660 


5.35 X 10"^ 


2.54 X 10 


-9 


4.46 X 10^ 


238 TT 
92 


0+ 


2+ 


E2 


44.910 


12.09 X 10"^ 


6.72 X 10" 


-10 


5.50 X 10^ 




5/2- 


7/2- 


Ml 


78.647 


3.07 X 10-3 


5.94 X 10 


-8 


1.81 X 10^ 


63 EU 


5/2+ 


7/2+ 


Ml 


21.532 


2.17 X 10"^ 


3.14 X 10" 


-10 


2.62 X 10^ 




7/2- 


9/2- 


Ml 


94.700 


6.76 X 10-3 


2.26 X 10 


-7 


7.42 X 10^ 



3 X 10-^° photons per nucleus per pulse for the case of the El transition of the 62^Sm isotope, two 
times larger than the one for the M 1 transition of 67^Ho. 

Generally, due to the very narrow widths of excited nuclear levels lying below 12.4 keV, the 
number of photons in the laser pulse that fulfill the resonance condition for the considered isotopes 
is small. The signal photon rates are for these cases up to ten orders of magnitude smaller than the 
ones for transitions above 12.4 keV. With regard to a possible experiment aimed at the observation 
of the direct laser-nucleus interaction, isotopes with low-lying first excited states below 12.4 keV 
have nevertheless the advantage that the laser can directly provide photons at the resonance en- 
ergy. Solid targets which have a very high density of nuclei can therefore be used for such an 
experiment. In this case, however, the presence of atomic electrons opens the IC decay channel 
for the nuclear excited state. The branching ratio for the nuclear decay can be related to the IC co- 
efficient a, defined as the ratio between the IC and radiative decay rates. Due to the small energies 



of the nuclear excited states, the IC decay is more probable, with a values between 



27 



28 



7.09 for the 



30 



9.396 keV transition of ^Kr and 920 for the 4.821 keV transition of ^^^Sm tm 
The photon signal per nucleus per pulse is then given by the fraction of the excited state population 
that will decay radiatively, Pee/{1 + ct). In order to estimate the signal photon rate per second, we 
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FIG. 1 : (Color online) Number of signal photons per nucleus per laser pulse A'^signal for several isotopes 
with first excited states below 12.4 keV (green squares) and above 12.4 keV (red crosses). The results are 
plotted versus the atomic number Z. The European XFEL considered has a pulse duration of 100 fs and an 
average brilliance of 1.6x 10^^ photons/(sec-mrad^-mm^-0.1%bandwidth). 



assume a focal diameter of the laser of 20 fxm together with a repetition rate of 40 kHz [[si |9[| and 
a density of 10^° nuclei/cm^ in the solid state target. In Figure l2t A) we present the photon signal 
per second for the M 1 and El nuclear transitions presented in Table HI The largest signal rate of 
80.75 photons/second corresponds to the 9.746 keV Ml transition of jf^Os. 

It should be noted, however, that the presence of the electrons in the solid target complicates 
the identification of laser driving of nuclei. The photon scattering off of the crystal electrons 
creates a considerable background. Furthermore, the irradiation of a strong XFEL laser on a 
solid target would lead to the formation of a plasma, invoking many other nuclear excitation or 
decay mechanisms involving high-energy charged particles, such as nuclear excitation by electron 



capture [|39L 4^|41|] or nuclear excitation by electron transition [|42l1 . Nuclear excitation and decay 



mechanisms in plasmas have been the subject of several studies 11431 



45n which show the 



importance of taking into account the electron-nucleus coupling. Given the narrow widths of the 



low- 



ying nuclear excited states, solid state effects like inhomogeneous broadening of the nuclear 



line 1461] may play an important role. We would like to stress, however, that similar difficulties have 



been experimentally overcome in the case of Mossbauer spectroscopy using synchrotron radiation 



[LZI], which led to the observation of interesting effects such as nuclear superradiance [|47 
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FIG. 2: (Color online) Signal photons per second for laser driving of nuclear transitions with energies (A) 
below 12.4 keV and (B) above 12.4 keV. The energy of the nuclear transition is given on the abscissa. The 
green squares denote El transitions, the blue circles E2 transitions and the red crosses Ml transitions. The 
considered experimental laser and target parameters are discussed in Sec. lIII B] 



nuclear exciton echos [14911 



The driving of nuclear levels with energies larger than 12.4 keV requires acceleration of the 

target nuclei in order to match the photon and nuclear transition energy. This restricts the target to 

ion beams, which have a small density in comparison with solid targets. With 2.5 x 10^° particles 

in a bunch length of r = 50 ns in an ion beam of 2 mm diameter as target [50], the particle density 

yields 5.3x10*^ cm~^. For the new Synchrotron SIS 100 that will be built in the future at FAIR [51 1 

O 

the beam parameters yield a particle density of 10^^ ions/cm^ [!5O0- Using this value to estimate the 
number of signal photons per second we obtain 2.35 x 10'^ s~^ for the El transition in gl'^Sm and 
1.16 X 10^ s~^ for the Ml transition in ey^Ho. Here, the calculation of the photon rate per second 
assumes the matching of the ion and laser pulse repetition rates. The signal photon rates per second 
for the El, E2 and Ml transitions given in Table HI] are presented in Figure [2tB). In addition to 
the larger signal photon rates, the experimental background is also more advantageous for such 
transitions, since fewer competing effects are present. In the case of bare ions, for instance, the IC 
nuclear decay channel can be suppressed completely. 

For the case of laser driving of the transition to the first excited state of gy^Ho, the boosted 
XFEL intensity of approximately X = 4 x 10^'' W/cm^ corresponds to an effective intensity Ics = 
10^^ W/cm^. As a result of laser driving of the nuclear transition, signal photons corresponding to 
the radiative decay of the nuclei can be observed. The dynamics of the nuclear few-level system 
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strongly depends on the laser intensity. For the considered XFEL laser intensity, the nuclear 
system remains on average almost in the ground state (the excitation rate per nucleus per laser 
pulse is Pee = 1-5 X 10~^° for gy^Ho), such that only few nuclei are excited in each laser pulse. 
With increasing intensity the nuclear excited state population starts to oscillate rapidly at about 
X = 10^^ W/cm^. A TT pulse that would directly transfer all nuclei to the excited state with 
no further oscillations can be found around X^, = 10^'' W/cm^. In analogy to atomic quantum 
optics, at X = 10^^ W/cm^, when the Rabi frequency equals the transition frequency, strong-field 
effects beyond the rotating wave approximation are expected to occur. However, note that with 
driving fields approaching the Schwinger intensity of about 10^^ W/cm^, other processes such as 
pair creation and vacuum screening become relevant and may prevent observation of the nuclear 
excitation. 

Finally, we would like to address as an interesting case study the possible excitation of the 
isomeric state of go^Th that lies only 7.6 eV above the ground state ll52ll . The metastable state is 
assumed to decay to the ground state by an Ml transition with a lifetime of approximately 5 h. 
While this lifetime corresponds to a very narrow width of the nuclear excited state, one may expect 



that in solid targets homogeneous line broadening can be of the order of 10 eV 



tunately, since the decay of the isomeric state has not been observed directly [|55 



3 



57 



151154]. Unfor- 



M, the 



reduced transition probability is unknown, thus making it impossible to give a realistic estimation 
of the laser-nucleus interaction matrix element. Detection of the laser-induced transition can be 
achieved by probing the hyperfine structure of a transition in the electronic shell, as it has been 
proposed in Ref. [|59|] . Since the energy value of the 229mrj]^ isomer was recently corrected from 



3.5 eV to 7.6 eV, numerous studies considering laser-assisted or 



nucleus and the electronic shell (see, for instance, Refs. [60, 



aser-driven coupling between the 



62L 



631) need revision. 



IV. SUMMARY 



We have considered the direct interaction between strong laser fields and nuclei focusing in 
particular on the driving of electric dipole-forbidden transitions. The main motivation for investi- 
gating such transitions is the limited variety of stable nuclear systems with El transitions in the 
interesting low-energy range. While El transitions with first excited states lying below 100 keV 
exist, the systems for which the ground state is stable or long-lived are limited to only a few cases. 
We have investigated laser-driven Ml and E2 transitions for a number of stable systems with 
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energies above and below the maximum photon energy of 12.4 keV envisaged at the European 
XFEL, and have compared the results with corresponding El cases. We found that on average, 
laser-driving of Ml transitions is as effective as the driving of -El transitions. In particular for 
excited levels lying below 12.4 keV, Ml transitions are the most promising for laser excitation. 
This is in contrast to atomic quantum optics, where higher multipole transitions are strongly sup- 
pressed. This equivalence of different multipole transitions in direct nucleus-laser interactions 
considerably enhances the range of suitable isotopes. Our analysis further shows that generally, 
in the search for a suitable experimental candidate system, the strength of the nuclear transition, 
usually given as the reduced transition probability, is a key quantity in estimating the magnitude 
of the effect. Finally, we discussed two possible scenarios for the experimental observation of 
direct laser-nucleus interaction. For transition energies below 12.4 keV, no target acceleration is 
required, such that high-density solid targets are preferable. We found signal rates of up to 80 
photons/second for the considered isotopes with low transition frequencies. Plasma effects oc- 
curring in solid state targets are, however, likely to complicate the experimental observation and 
require an adequate theoretical description. Above the maximum source photon energy, the use of 
accelerated ion beams is required. These have a lower target density than solid-state targets, but 
also a much reduced background. Using realistic parameters, we found for this case signal photon 
rates of up to lO'' photons per second. 
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